Activated dynamics and effective temperature in a steady state sheared glass 
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We conduct nonequilibrium molecular dynamics simulations to measure the shear stress a, the 
average inherent structure energy (£rs), and the effective temperature T c g of a sheared model glass 
as a function of bath temperature T and shear strain rate 7. For T above the glass transition 
temperature To, the rheology approaches a Newtonian limit and T c ff — > T as 7 — * 0, while for 
T < To, a approaches a yield stress and T e ff approaches a limiting value near To. In the shear- 
dominated regime at high T, high 7 or at low T, we find that a and (-Bis) each collapse onto a single 
curve as a function of T e «. This indicates that T c s is controlling behavior in this regime. 

PACS numbers: 05.70.Ln, 64.70.Pf, 83.50.Ax 



When a liquid is quenched through the glass transition 
temperature To, it falls out of equilibrium and becomes a 
glass. Below To, the thermal energy is insufficient for the 
system to surmount energy barriers on accessible time 
scales, so the glass explores relatively few configurations. 
However, if the glass is held in contact with a thermal 
reservoir and sheared at a fixed rate, it can reach a steady 
state in which it explores many different minima in the 
energy landscape, even if the temperature of the reservoir 
is well below To [l| . The steadily sheared glass is far from 
equilibrium: energy is continually supplied on long time 
and length scales via the boundaries, and is removed on 
short scales by the thermal reservoir. Nonetheless, sim- 
ulations show that fluctuations in such systems are well 
described by an effective temperature T e g that is higher 
than the bath temperature 0, H, H|, as predicted theo- 
retically [f|. Nine different definitions yield a common 
value of T e fi H, 0, @, 0, @| , providing strong numerical 
evidence for the utility of the concept. Most recently, Ilg 
and Barrat 4] showed that T c g controls activated transi- 
tion rates of a test probe consisting of a dimer connected 
by a double- well potential, embedded in a sheared model 
glass. The rate of crossing the energy barrier separat- 
ing the two wells has an Arrhenius form cxp(— AE/T c g), 
where T e g is consistent with previous definitions. 

In this paper, we shift the focus from testing the valid- 
ity of effective temperature to examining its importance 
for materials properties. We show that T e g plays a crit- 
ical role in fluidizing a glass. We find that the shear 
stress collapses onto a single curve depending only on 
T e g whenever the shear rate is high enough to dominate 
over thermal effects. The average inherent structure en- 
ergy collapses in similar fashion. These findings suggest 
that T c ff activates particle rearrangements necessary for 
flow, much as thermal fluctuations do in an equilibrium 
liquid, supporting the idea that a common framework 
might describe unjamming by mechanical forcing and by 
temperature 0. 

Our simulation model is a two-dimensional glass- 
forming liquid composed of disks interacting via a purely 
repulsive harmonic potential [3j|. We study 50:50 mix- 



tures of disks of diameter ratio 1 : 1.4 and equal mass. 
The area fraction is fixed at <\> = 0.9. Most of the results 
are based on simulations of 400 disks, but we carried 
out simulations of up to 6400 disks to confirm that none 
of these results have any appreciable system size depen- 
dence. Units in this paper are measured with the smaller 
particle diameter, the interaction spring constant, the 
particle mass, and the Boltzmann factor set equal to 1. 
This yields a unit time period on the order of a binary 
collision time. We uniformly and steadily shear the sys- 
tem at a strain rate 7 and couple the system to a heat 
bath at temperature T by solving the SUod equations 
of motion with Lees-Edwards periodic boundary condi- 
tions and a Gaussian thermostat ■ We integrate these 
equations using a fourth-order Gear corrector-predictor 
algorithm with a time step of 0.01. We obtain the same 
results with a Nose-Hoover thermostat. We use between 
five and twenty simulation runs for each set of (T, 7) . For 
each simulation, we collect data over at least nine strain 
units after an equilibration period of several strain units. 

We measure T e fr from the relation between the static 
linear response and the variance of the pressure 0, [HI • 
In equilibrium at fixed N, T, and A, this relation is 



(i) 



where p is the pressure, A is the area, and x is the hy- 
pervirial as defined in [l0(. T e g is defined by replacing 
T with T c ft in the left hand side of Eq. |T]). We mea- 
sure d(p)/dA by running simulations at (f> = 0.897 and 
4> = 0.903, using the same protocol and a similar quan- 
tity of simulations as for <f> = 0.9. Measurements from 
this definition have been compared to those from many 
other definitions of T e g for a zero-temperature sheared 
foam [jj • We find that T e fj is consistent with less precise 
measurements of T e w defined by the Green-Kubo relation 
for shear viscosity 0, O'Hern et al. 0] showed that 
T e ff from pressure fluctuations agrees with that derived 
from the time-dependent linear response of density fluc- 
tuations 0] over a range of parameters for the system we 
use. We also find consistency between these two defini- 
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FIG. 1: a (a) and T e s (b) vs. 7 for several values of T and 7. 
In all figures, dashed lines and open points indicate T > To. 
Solid lines and solid points indicate T < To. In (a), the 
straight line indicates a slope of 1, expected for a Newtonian 
fluid. In (b), the horizontal line indicates To, while the iso- 
lated points near the left axis indicate the values of T. 



tions, except at very low strain rates where the logarith- 
mic time-dependence of the diffusivity, expected in two 
dimensions [11] , is apparent within the time scale of the 
density correlation function. This long-time tail does not 
affect the viscosity because at the high area fraction and 
low temperatures studied, the kinetic contribution to the 
viscosity is much smaller than the potential contribution. 

Figure [1] shows the measured quantities, stress a and 
effective temperature T e g, as functions of bath tempera- 
ture T and strain rate 7. Figure [lj a) shows a bifurcation 
of a about a critical temperature Tq = 0.0012 at low 7, 



similar to what has been observed in experiments 12 1 
and earlier simulations 0] of sheared glass. Bifurcations 
about a critical volume fraction have been observed in 
experiments on colloids [13] and emulsions fl4ll and in 
simulations of frictionless spherical packings [151 ]. At To, 
the stress obeys power-law scaling with 7, with a ~ 7 3 . 

The two sides of the bifurcation in Fig. HJa) describe 
two different low-7 limits. On the T > Tq side, er ap- 
proaches a ~ 7 at low 7, suggesting that the shear vis- 
cosity 77 = cr/7 enters a Newtonian regime for T > Tq and 



sufficiently low 7. For T > 0.0017, we have reached strain 
rates low enough that 77 becomes independent of 7. We 
define the equilibrium viscosity 7y eq as the shear viscos- 
ity in this Newtonian regime. For T < T < 0.0017, 
we do not reach the Newtonian regime at accessible 
time scales/strain rates. However, we find that for all 
T > 0.0015, 77(7) is well-described by the phenomenolog- 
ical Ellis equation [B], l/r]{a) = l/r) + TO -i/™ CT (i-«)/™ ; 
which interpolates between Newtonian and power law 
rheology and allows us to define r] e<1 down to T = 0.0015. 

In contrast, for T < Tq we observe apparent yield stress 
rheology on the time scale of our simulations [l7j. We 
find a — Cyioid oc 7^ over the lowest two decades of 7 that 
our simulation can access, with the exponent b ranging 
between 0.4 (T = 0.0012) and 0.6 (T = 0.0001), similar 
to that observed in emulsions G3 • This implies that the 

eld- 



viscosity diverges as 77 = a /(a — tTyicid) 1 ^ as a — > a 



Figure QJb) shows the dependence of T c g on T and 7. 
For all T, T e s approaches a limiting value in the qua- 
sistatic limit, 7^0. For T > Tq, that limiting value is 
simply T. However, for T < To, T e g appears to saturate 
to a value T e g t Q near To Q. Such a saturation of T c g at 
low 7 has been observed in experiments on sheared gran- 
ular packings 18, [l^]. This apparent quasistatic limit 



suggests that T c s,o is a property of the unsheared glass 
at bath temperature T, describing the disorder associated 
with different minima in the energy landscape [20j | . 

In Fig. [2]^a), we compare the approach to jamming as 
T — > To with the approach at fixed T < Tq and 7 — > 0, 
parameterizing the latter approach by T e g rather than 7. 
The dependence of 77 on 1 /T e ff is similar to that of r] eq on 
1 /T [21| , but we find no collapse among the different ap- 
proaches to jamming. Along the equilibrium approach, 
77 has the Arrhenius form 77 — 7700 exp^^/T) at high 
T. The non-equilibrium approaches exhibit no Arrhe- 
nius regime in 77. However, another reasonable measure 
of relaxation time, r s hear = 7 (not shown), does vary 
in Arrhenius fashion with T e ff at high T e g . Along all ap- 
proaches, 77 is super- Arrhenius for T c s near To or T e gfl. 

Figure f^b) demonstrates that the stress a collapses as 
a function of T c g in the shear-dominated regime. We find 
that <7 varies between two limits, depending on T and 7. 
For T > Tq and 7^0, the shear stress approaches the 
Newtonian limit, a = 77,^7, while T e ff approaches T. In 
this regime, thermal fluctuations T dominate over shear- 
induced ones. However, for T <c To and any 7, or for 
T>To and high 7, the system crosses over to another 
regime, where shear-induced fluctuations dominate over 
thermal ones. In this shear-driven regime, a depends on 
T and 7 only through T e g. The dependence of a on T e g 
follows the simple form 



a w cr exp(-A J B/T eff ), 



(2) 



represented by a dotted line in Fig. (2^b). The col- 
lapse of the data from different T and 7 onto this same 
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FIG. 2: (a) f] vs. 1/Teff for several values of T and 7. Asterisks 
correspond to the Newtonian viscosities (7 — > 0, T e g = T). 
(b) a vs. 1/T e a. The straight dotted line is a fit to the form 
a — aoexp( — AE/T e ff) for the three lowest bath tempera- 
tures, (c) (-Bis) vs. 1/T e g. Asterisks indicate (Eis) for 7 = 
(T«ff = T). 



curve indicates that the dynamics are most properly de- 
scribed as being controlled by T c ff, not 7 or T. In 
the shear-dominated regime, the height of the energy 
scale, AE = 0.0028 ± 0.0001, is independent of T and 
7. Its value is comparable to the activation energy 
E A = 0.0023 ± 0.0002 extracted from the high-T equi- 
librium Arrhenius relationship 77 = rjao bxp(Ea/T). The 
correspondence between a and T c g suggests that the ex- 
istence of a finite yield stress would imply a nonzero qua- 
sistatic value of T e s t Q. 

Equation ^ suggests a simple scenario. Suppose that 
the power per unit area supplied globally by shearing, 
cry, were used to overcome local energy barriers of height 
AE at a rate of R transitions per unit area per unit time. 



Using cry = RAE and the observed Eq. ((2|), we find 

R = R jcxp(-AE/T cS ), 



(3) 



with i?o = ao/AE = 4 ± 1. This implies that the rate 
of barrier crossing is Arrhenius in T g with an attempt 
frequency proportional to the strain rate 7. Local re- 
arrangements whose rate scales with 7 are observed in 
foam experiments [22| and simulations [23| , as well as in 
athermal quasistatic simulations of soft disks 



24] 



Although it was derived from a steady-state relation- 
ship, Eq. [3] motivates a hypothesis for the time evolution 
of T e ff. Such equations [5^, [26| have been used in the 



context of shear transformation zone (STZ) theory [27 1. 
Following previous work [25|, we suppose that the equa- 
tion should balance a heating rate proportional to the 
work done on the system with a relaxation rate propor- 
tional to exp(—Ei/T c g). As in Eq. [3]and differing from 
[25| . we suggest that the relaxation rate should depend 
on a scalar rate at which the system explores new con- 
figurations. These considerations yield 



T eff oc Q - vexp(-E 1 /T efl ), 



(4) 



where Q oc cy [28[ and the attempt frequency v oc 7 
in the case of steady-state shear. Equation [4] is consis- 
tent with Eq. [2] in the steady-state limit. For T e g near 
its steady-state value T e ff,sS) Eq. 0] reduces to T e g oc 
CT7(7rff,ss - Teff), as used in STZ theory 

The scenario suggested by Eqs. OH] is that T e g acti- 
vates the system over barriers whose height is indepen- 
dent of T ff. However, measurements of the average in- 
herent structure energy (Eis) of the system suggest 
that the height of energy barriers does depend on T e g. 
We measure Eis of the sheared system by taking config- 
urations explored during steady-state shear and quench- 
ing them to their local energy minima by the conjugate- 
gradient technique. For comparison, we also measure 
(-Eis) (T c ff = T) for the equilibrium system at 7 = 
above the glass transition temperature. Fig. shows 
that (Eis) (Teff = T) is flat at high T. For T < T MC , 
where Tmc marks the onset of super- Arrhenius behavior 
in Fig.[2fa), we find that {Eis)(T e g = T) decreases mono- 
tonically with decreasing T, in agreement with [30I ]. At 
low T, the system visits deeper potential energy basins, 
presumably separated by higher barriers, consistent with 
the super- Arrhenius viscosity 77 = rj 00 exp(EA(T)/T), 
where Ea(T) is the T-dependent barrier height. 

Figure [^c) shows that (Ei$)(T e g) of the sheared sys- 
tem is distinct from the equilibrium curve (2oT | but also 
decreases as T e g — > T e g t o. Moreover, the data appear to 
collapse in the shear-dominated regime, as in Fig. [2Kb) . 
At sufficiently low T and/or high 7, the data collapse 
onto a single curve, but cross over to the equilibrium 
curve at high T and low 7. For each T, the point at 
which (-Eis) (T e ff) begins to decrease corresponds to the 
upturn of viscosity, suggesting that energy barriers in- 
crease with decreasing (Eis). 
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In summary, there are two possible explanations for 
the upturn of r\ as T e g — > T c ff,o- The first explanation 
is that relaxation rates are controlled by T c ff-activated 
transitions over barriers whose heights increase with de- 
creasing T c ff. This view is supported by an STZ analysis 
of our data 3l| . The second explanation is that the bar- 
rier heights overcome during the shearing process do not 
depend strongly on T e fj, even though the energy minima 
decrease with decreasing T c g. In that case, the super- 
Arrhcnius behavior is due to Eq. [2] and the divergence 
in r\ w (Tyicid/7 as 7 — > 0. Further study, particularly of 
transients, is needed to resolve this issue. 

Finally, we revisit the issue of the validity of the ef- 
fective temperature concept. Nine independent defini- 
tions of temperature have been shown to yield consis- 
tent values of T e g, within numerical error: the relation 
of density @, [1] and pressure fluctuations Q at nonzero 
wavevectors to the associated response at long times; the 
relation of fluctuations in the total pressure, stress, and 
energy to static response [3L [A\ the Einstein relation be- 
tween diffusion and drag [2, lOjLZI L the derivative of en- 
tropy with respect to energy [y, Ijj ; the fluctuations of 
a low- frequency harmonic oscillator Q; and the barrier 
crossing rate of a test two-level system jij. However, 
there are definitions that do not yield consistent values 
of T e g: the relation of fluctuations in the total devia- 
toric pressure and in the vorticity component of pressure 
to the response at long times [3(. Thus, the concept of 
effective temperature, even when restricted to long-time- 
scale properties, is only approximate [jl. l32|. 

If the concept of effective temperature is not rigorously 
valid for sheared glasses, why would it be of any interest 
for these systems? Our results provide an answer: the 
effective temperature critically affects materials proper- 
ties by setting the energy scale for fluctuations that kick 
flowing glasses over energy barriers. 
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